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Using a sample of 122 million T(3S') decays collected with the BABAR detector at the PEP-II asymmetric 
energy collider at the SLAC National Accelerator Laboratory, we measure the ratio Rt/j, = BR{T(IS) — ► 
T~^T~)/ BR{T{IS) — > the measurement is intended as a test of lepton universality and as a possible 

search for a light pseudoscalar Higgs boson in Next to Minimal Supersymmetric Standard Model (NMSSM) 
scenarios. Such a boson could appear in a deviation of the ratio Rt/j, from the Standard Model expectation, 
that is 1, except for small lepton mass corrections. The analysis exploits the decays T(35) — + T(lS')7r+7r~ , 
T(1S) ^ i+r , where / = /i, r. 



1. Introduction 

In the Standard Model (SM), the coupHngs between 
gauge bosons and leptons are independent of lepton 
flavor. Therefore, the branching ratio BR{T{1S) — > 
should basically not depend on the lepton I con- 
sidered. The quantity: 

^ BR{T{IS) ^ l+n 

BR{T{1S) ^ I'+l'-)' ^ ' 

where 1,1' = e,^,T, is thus expected to be very near 
to one, the discrepancy from unity being due to small 
lepton-mass effects (for instance, the greatest correc- 
tion should turn out in R^^ ~ 0.992). 

Beyond the SM, deviations of Rw from this expec- 
tation are possible. In particular, in the Next to Min- 
imal Supersymmetric Standard Model (NMSSM) [T] 
scenario, there is the hypothesis of existence of a light 
pseudo-scalar Higgs boson, namely A^, which could be 
escaped to LEP bounds [2J [3] . For instance, consid- 
ering the leptonic decays of the T(15), such a boson 
could mediate the decay chain of the resonance in one 
of the following modes [H |5] : 

T{1S)^ A°-f,A° ^l+r (2) 

or 

T(l^) ^ 7jb{lSh,7jb{lS) ^° i+r. (3) 

The latter decay implies a mixing between r]b{lS), i.e. 
the ground state of the bottomonium family, recently 
discovered by BaBAR [6] , and which does not have the 
right quantum numbers to decay in a pair of leptons, 
and the light Higgs A'^ . If the photon was present but 
remained undetected because of its softness, the lep- 
ton pair would be ascribed to the T(IS'). This would 
result in a deviation of i?;// from the SM expectation, 
that is in a lepton universality violation. Similar con- 
siderations hold on for the higher-mass T resonances 
too. The coupling between A*^ and the lepton being 
proportional to the lepton mass, this effect should be 



more evident when one of the leptons is a r. Theoret- 
ical calculations [U [S] foresee effects up to an order of 
- 10%. 

The CLEO Collaboration has studied the decays 
T{nS) /^^A* > T+T^, with n = 1,2,3, using a sam- 
ple of about 1 fb^^ of data taken at each T{nS) peak, 
corresponding to about 10^ T resonances. In particu- 
lar, for the T(IS'), they measured the branching frac- 
tion ratio: 

Rr^, = 1.02 ± 0.02(stat.)±0.05(syst.), 

deriving: 

B{T{1S) ^ ?767) X B(% ^ A" ^ T+T-) < 0.27% 
at the 95% of confidence level [7j. 

2. The BaBar detector and the data 
sample 

We present a measurement of Rt^, in T(IS') decays 
exploiting the sample of T(3S') decays collected by the 
BaBar detector at the PEP-II asymmetric e+e" col- 
lider at the SLAC National Accelerator Laboratory. 
The available statistics (^ 28 fb^^, corresponding to 
(121.8 ± 1.2) X 10^ events) represents the largest sam- 
ple of T(35) decays ever collected. The 2.4 fb~^ of 
data collected 30 MeV below the T(3S') resonance 
(off-resonance sample) are also used for background 
studies. 

The BaBar detector is described in detail else- 
where |H]. Charged-particle momenta are measured 
in a tracking system consisting of a five-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer cen- 
tral drift chamber (DCH), both situated in a 1.5-T 
axial magnetic field. Charged-particle identification 
is based on the dE/dx measured in the SVT and 
DCH, and on a measurement of the photons pro- 
duced in the synthetic fused-silica bars of the ring- 
imaging Cherenkov detector. A CsI(Tl) electromag- 
netic calorimeter is used to detect and identify pho- 
tons and electrons, while muons are identified in the 
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instrumented flux return of the magnet. We use the 
GEANT [9] software to simulate interactions of par- 
ticles traversing the BaBAR detector, taking into ac- 
count the varying detector conditions and beam back- 
grounds. 



3. Analysis strategy 

The signal events (T(IS') /^^/-^^ or 

r+T~) are tagged through the transition 
T{3S) T(lS')7r+7r^, whose branching fraction 
is BR{T{3S) -> T(lS')7r+7r-) - 4.8% [TD]. We 
consider all the decays of r to one charged track plus 
undetected neutrals: this choice leads to a very clear 
environment for our study, since both the typologies 
of final states (/i^/i^ and r+r") contain four charged 
tracks, allowing for a better handling of backgrounds 
and a good trigger performance. 

Because of the presence of neutrinos in the final 
state, the T(35') ^ T(lS')7r+7r-, T(IS') ^ t+t- de- 
cays are not fully reconstructed, and suffer from a 
larger contamination from background with respect to 
the T{1S) M^M^ events. For this reason the event 
selections are mostly different for the two categories 
of final states. 

The event selections are tuned using Monte Carlo 
simulated samples. The main sources of background 
are e+e^ qq (with q = u,d, s, c) and e'^e~ t^t^ 
continuum events, Bhabha events and generic T(IS') 
decays; the latter produce a peaking yield. The two 
typologies of signal events {^'^fj,~ and r+r^ final 
states) are separated thanks to a cut on the differ- 
ence between the visible and the reconstructed energy, 
calculated in the e+e~ center of mass frame (CM): 
indeed, no missing energy is expected in T(IS') 

events, that are completely reconstructed, while 
an amount of missing energy is actually present in 
T(IS') T^T~ decays. 

Moreover, to further reduce the background in the 
T+r~ final states, a boosted decision tree [TT] exploit- 
ing several shape and kinematic variables is used. The 
performance of the classifier is assessed using signal 
Monte Carlo simulations and off-resonance data. Af- 
ter the application of the multivariate approach, the 
separation between signal and background is sensi- 
tively improved, as visible in Figure [l] 

Additional cuts are applied for the /x^/i^ selection 
on the cosine of the angle between the two-lepton can- 
didates calculated with respect to the CM frame, for 
the r+r~ selection on the difference in the energy of 
the T(3S') and the T(IS'), and for both the selections 
on the momentum of each of the two pions, as well as 
of the pair of them. 

The signal extraction efficiency is e^^ ~ 45% {err ~ 
17%) for n'^iJ,^ (t+t^) final states. 



1 MVA output (or method: BDtI * TMVA 




BDT 

Figure 1: Output variable resulting from the multivariate 
analysis method, exploiting a boosted decision tree as a 
discriminator tool, for signal Monte Carlo simulation (solid 
blue) and for off-resonance data (striped red). 

4. Signal extraction 

The signal extraction is performed using an ex- 
tended and unbinned maximum-likelihood fit; the 
likelihood for the T(IS') f^~^IJ-~ events is 2- 

dimensional and based on the variables AM = 
Mr(3S) — Mx{is), that is the difference between the 
invariant masses of the two T resonances involved in 
the decay chain, and M^+^- , the invariant dimuon 
mass. Instead, for the T(IS') t~^t~ events, the sig- 
nal extraction cannot be performed in the same way: 
the partial reconstruction of the final state, due to the 
presence of neutrinos from the r's decays, means that 
neither AM nor the dilepton invariant mass can be 
used to discriminate between signal and background. 
Therefore, we extract the signal yield performing a 1- 
dimensional fit to the recoil mass of the pair of pions 
involved in the decay chain, defined as: 

m:%':- - + - 2v:^./mi~^pi^^,„ (4) 

where y/s indicates the nominal CM energy [10], Mtttt 
indicates the invariant mass of the two pions, and 
PttttCM is the momentum of the two pions calculated 
in the CM frame. Because this variable is related only 
to the T(3S') 7r+7r^T(15) transition, it cannot dis- 
tinguish between T(15') —^ t~^t~ events and other 
T(IS') decays or the Higgs-mediated events of Eq. [2] 
and Eq. [S] We use a simulated sample of inclusive 
T(IS') decays to verify that this background becomes 
negligible when the event selection is applied. Ad- 
ditional tests are performed with the data sample to 
prove this conclusion, as discussed in detail below. An 
enhancement of the T(IS') — > r+r^ yield would thus 
constitute a possible signal of New Physics. 
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The two fits are performed simultaneously, and the 
value of Rrfi is returned. The likelihood is written as 

Ce.t^ 11^* (5) 

i=l 

where Vi is: 

+ '—N^Rr^vr{M:T:-) + 

+ A^,fe,.7'f'^(K?;-) (6) 

and Nfj_, Ni,kg^ and Nfj^g^- indicate the number of the 
signal events extracted in the fJ-^ fJ-" sample, the num- 
ber of background events extracted in the fi~^ n~ sam- 
ple and the number of background events extracted in 
the T~^T~ sample, respectively. 

In the fit, the probability density functions (PDFs) 
for the signal components in T(15) /^^M^ sample 
are taken from a subset of data (about one tenth of 
the complete available statistics, which is discarded 
from the eventual result in order to avoid any possi- 
ble bias); the AM distribution is described by a triple 
Gaussian function with uncommon means, while the 
M^+^i- distribution is modeled by an analytical func- 
tion approximating a Gaussian with mean value equal 
to ^ and different left and right widths a{L,R), plus 
asymmetric non gaussian tails a(L,R), defined as: 



= e.p{ - 2^,^^^ ^"^^^^ ^^^^ _ J. (7) 

In T(IS') T^T" events, the PDF for the signal 
component is taken from the data sample T(IS') — s- 
sample (where the dipion recoil mass is not used 
as a variable for the signal extraction and is therefore 
independent of the fit procedure); the M^!j:™_ distri- 
bution is described by the function defined in Eq. [7) 
All the parameters of the signal PDFs are fixed in 
the fit, while the background shapes (flat functions 
in /x"*"/!" sample and polynomial function of first or- 
der in T+r^ sample) are fioated in the fit. We use 
the off-resonance sample to verify that the analytic 
description of the background shapes is appropriate. 

The result of the simultaneous fit procedure is 
shown in Figure [2] while Table |l] shows the results 
for the signal and background yields, along with the 
results of the jd.o.f. tests, made to compare the 
data distribution with the curve chosen to describe it 
(where d.o.f. indicates the number of degrees of free- 
dom, namely equal to the number of histogram bins 
minus the number of parameters fixed in the PDF). 



Table I Results of the simultaneous fit: signal and back- 
ground yields; -Rt^j; output of /d.o.f . tests. The results 
are only statistical; no systematic correction is applied yet. 





60192 ± 269 




3136 ±92 




44023 ± 253 




1.006 ±0.010 


X^{AM)/d.o.f. 


176/92 


X\M^+^-)/ d.o.f. 


133/95 


x\M;r°^)/d.o.f. 


158/95 



5. Systematic uncertainties 

In the measurement of Rr^,, several systematic un- 
certainties are common to the numerator and denom- 
inator of the ratio, and therefore cancel out. These 
include uncertainties on the luminosity, the T(3S') 
production cross section, the T(35) T(15)7r+7r^ 
branching fractions, as well as uncertainties from dif- 
ferences between data and simulation in track recon- 
struction efficiencies and the common event selection. 

The residual systematic uncertainties that do not 
cancel arise from: 

• systematic uncertainty related to the event se- 
lection efficiency; 

• differences between data and simulation in the 
particle identification (PID) efficiency; 

• differences between data and simulation in trig- 
ger and background filter efficiency; 

• imperfect knowledge of the signal and back- 
ground shapes used in the fit to extract the num- 
ber of signal events on data. 

The systematic uncertainty due to the event selec- 
tion is evaluated as the discrepancy between data and 
simulation in the variation of the efficiency produced 
by removing one cut at a time from the complete event 
selection. This procedure is performed for those selec- 
tion criteria which are not common between the two 
channels (and hence do not cancel) and also for the se- 
lection criteria which use different values of the same 
variables (which partially cancel). The systematic er- 
ror is 1.5%, mainly due to the effect of the cut applied 
to the output variable of the multivariate technique 
used for T(15) t^t~ event selection. 

The systematic uncertainty related to the difference 
between data and simulation in the PID efficiency only 
applies to T(IS') 1^'^ ^^~ events, and is estimated 
by exploiting two independent samples: one where 
both leptons are required to be identified as /it leptons, 
and another where exactly one final charged track is 
a muon. We determine the relative inefficiency for 
requiring the second sample instead of the first one. 
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Figure 2: Results of the fit procedure for the AA/ distribution (top left), for the M^+^- distribution (top right) and for 
the MJJ+'^°_ distribution (bottom). In each plot the dashed green line represents the background shape, while the solid 
blue line is the sum of signal and background contributions. The residuals of the fit are shown under each plot, with the 
solid red lines indicating the 2<j level of agreement. 



both on data and on simulation; the ratio of the two 
results gives an efficiency correction of 1.023 and a 
related systematic error of 0.6%. 

The systematic uncertainties due to the differences 
between data and simulation in trigger and back- 
ground digital filters (BGF) efficiency are small both 
in T(IS') fJ-~^ fJ-~ and in T(IS') t'^t~ events, and 
they cancel at least in part in the ratio. These values 
are separately estimated for each channel by compar- 
ing the efficiency for the signal extraction with and 
without the trigger requirements (and similarly with 
or without the BGF), both on data and on simula- 
tion. Concerning the trigger, a correction of 1.020 is 
needed for the e-^r efficiency, together with a system- 
atic uncertainty of 0.10% for T(15) — > T+r~ events, 
while a systematic uncertainty of 0.18% is quoted for 
T(IS') — *■ events. The contribution given by 

BGF efficiency to the total uncertainty is negligible. 

The imperfect knowledge of the signal and back- 
ground shapes used in the signal extraction procedure 
is taken into account considering: 

• the systematic uncertainty resulting from fixing 
the parameters of the functional forms describ- 
ing the signal distributions (up to 1.7%); this 



error is estimated by repeating the fit procedure 
changing of ±lcr all the parameters fixed, one 
at a time, and summing in quadrature the dis- 
crepancies between the nominal result (Table |l]) 
and the one obtained in each of the modified 
configurations; 

• the systematic uncertainty due to the choice of 
the background parameterization for each distri- 
bution used in the fit procedure: it is evaluated 
by using alternative parameterizations describ- 
ing the background shape (polynomial functions 
of first order for AM and M^+^~, flat func- 
tion and polynomial function of second order 
for M''+'^°_ ) and summing in quadrature the dis- 
crepancies between the nominal result and the 
one obtained in each of the modified configura- 
tions (up to 0.28%); 

• the systematic discrepancy between the sam- 
ple used to determine the parameters of each 
distribution. For T(IS') — > /.t^^~ events, a 
sub-sample of data themselves is used, lead- 
ing to no expected bias and therefore without 
a systematic assigned. For T(15) t^t^ 
events, the fJ-^fJ-'' sample is exploited, and an 
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estimate of the systematic discrepancy between 
the two samples is needed. It is evaluated by 
re-weighting the parameters for the M''+'^°_ dis- 
tribution with the ratio between the parameters 
obtained comparing the fit result on the Monte 
Carlo samples for both the leptonic T{1S) de- 
cay channels. We estimate this contribution to 
be about 0.10%. Moreover, in order to consider 
further discrepancies possibly not taken into ac- 
count by simulations, an additional systematic 
uncertainty is estimated by repeating the fit pro- 
cedure while allowing the global width of the 
function defined in Eq. |7] (which describes the 
M^%'^°_ shape for the signal) to float; the result 
is a discrepancy of 0.05% in Rt^. The total con- 
tribution is therefore 0.11%. 

Finally the statistics of the Monte Carlo simulated 
samples used gives a contribution to the systematic 
uncertainty of less than 0.1% in both the leptonic final 
states. 

Additional studies have been performed in order 
to check the possible background contribution to the 
T+r~ sample after applying the selection. The sig- 
nal extraction is repeated requiring the two r leptons 
in the T{3S) T(lS')7r+7r~ r+T~7r"'"7r~ decay 
chain to be detected in the final state e+/i^ or /i+e~. 
The central values for Rr^ obtained with the two pro- 
cedures on the control sub-sample of the complete 
dataset (already used in other steps of the analysis 
and discarded from the final result) are compatible: 
= 1.006 ± 0.029(stat.) and = 0.999 ± 

0.044(stat.), where the first value indicates the result 
obtained with the nominal procedure and the second 
one the value obtained requiring r+r^ e*/i^+ neu- 
trinos. This further confirms that the peaking back- 
ground due to generic T(IS') decays is indeed negligi- 
ble. 

The summary of the systematic uncertainties is 
shown in Table |TTj The total systematic error, ob- 
tained by summing in quadrature all the contribu- 
tions, is estimated to be 2.4%. 



6. Results 

Including all the systematic corrections, the prelim- 
inary result obtained by BaBAR is: 

Rr^ = 1.009 ± 0.010(stoi.) ± 0.024(sysi.). 

Therefore we do not observe any significant devia- 
tion of the ratio Rrfi from the SM expectation. With 
the present result the BaBAR Collaboration achieves a 
greater level of precision with respect to the previous 
best value [7], improving both the statistical and the 
systematic uncertainties. 



Table II Summary of systematic errors and correction 
to the efficiency for /i"*"^" channel (first column) and for 
t'^t~ channel (second column). The total systematic un- 
certainty is also shown. 



Systematic error: 




T+r- 


Event selection 


1.5% 


PID 


0.6% 




Trigger 


0.18% 


0.10% 


BGF 


negl. 


negl. 


PDF parameters 


1.7% 


Background PDFs 


0.28% 


Agreement vs. r+r^ in i\f^+'5^_ 




0.11% 


MC statistics 


0.08% 


0.09% 


TOTAL 


2.4% 



This result is still preliminary, since a further im- 
provement in the systematic precision is on going. The 
finalization of the analysis is expected for soon. 
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